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A nondissipative supercurrent state of a Josephson junction is metastable with respect to the 
formation of a finite-resistance state. This transition is driven by fluctuations, thermal at high 
temperatures and quantum at low temperatures. We evaluate the life time of such a state due to 
quantum fluctuations in the limit when the supercurrent is approaching the critical current. The 
decay probability is determined by the instanton action for the superconducting phase difference 
across the junction. At low temperatures, dynamics of the phase is massive and is determined by the 
effective capacitance, which is a sum of the geometric and intrinsic capacitance of the junction. We 
model the central part of the Josephson junction either by an arbitrary short mesoscopic conductor 
described by the set of its transmission coefficients, or by a diffusive wire of an arbitrary length. 

The intrinsic capacitance can generally be estimated as C* ~ G/Eg, where G is the normal-state 
conductance of the junction and Eg is the proximity minigap in its normal part. The obtained 
capacitance is sufficiently large to qualitatively explain hysteretic behavior of the current-voltage 
characteristic even in the absence of overheating. 

PACS numbers: 74.40.Gh, 74.45.-l-c 


I. INTRODUCTION 

Macroscopic quantum tunneling is a fascinating man¬ 
ifestation of quantum-mechanical behavior in large-scale 
systems with many degrees of freedom. It is responsible 
for a finite life time of a metastable state that cannot 
decay classically by thermal activation at zero temper¬ 
ature. Since 1980-ies, macroscopic quantum tunneling 
has been studied in a number of condensed-matter sys¬ 
tems: various types of Josephson junctionsE^' Joseph¬ 
son junction arrays,^ phase-slip centers,vo rtices 
in superconductors,^ small ferromagnetic particlesj ^^ * ^^ ! 
etc. 

Theoretical description of macroscopic quantum phe¬ 
nomena is based on the concept of a collective coordi¬ 
nate, X- In. the simplest cases, it can be considered as 
a slow variable, which allows to integrate out the other 
electronic degrees of freedom and end up with an effec¬ 
tive quantum mechanics for x(t). The resulting dynam¬ 
ics of the collective degree of freedom is generically non¬ 
local in time since interaction with other modes produces 
the retardation effect, intimately related with dissipation. 
Intensive studies of dissipative quantum mechanic^I^H^oi 
have been tigered by the pioneering work of Caldeira 
and Leggett .El 

Among various systems, the Josephson junction can 
be considered as a prototypical model of macroscopic 
quantum tunneling. Here the superconducting phase 
difference across the junction, x^ pl&ys the role of a 
collective coordinate. A nondissipative Josephson cur¬ 
rent can run through the system, described by a certain 
current-phase relation I{x)¥^ nonsinusoidal in general 
(Fig. [^), with the critical (maximal) current Ic reached 
at some phase difference Xc- The current can be ob- 


a) b) 



FIG. 1. (a) Typical current-phase relation for a Josephson 

junction, (b) Free energy F{x) of a current-biased junction 
vs. the phase difference Xi for the current / slightly below E- 


tained by differentiating the free energy of the junction: 
I{x) = {2e/h) dFQ{x)/dx- In the current-biased regime 
with the driving current I, the equilibrium states (x = Xi 
mod 27r) correspond to the minima of the Legendre- 
transformed free energy, F{x) = Fb(x) —(Ii/2e)/x, which 
has the standard form of a washboard potential shown 
in Fig. The supercurrent state with X = Xi is 

metastable and does decay (due to quantum or thermal 
fluctuations) into a resistive branch. Then the junction 
usually stays in the dissipative regime unless / is de¬ 
creased to a smaller retrapping current, resulting in a 
hysteretic current-voltage characteristic.!^ 

Early studie s of the supercurrent decay in Josephson 
j unction assumed the tunnel limit, when the su¬ 
perconducting terminals are coupled through an insulat¬ 
ing layer without its own electron dynamics (SIS junc¬ 
tions). Possible dissipation and charging effects were in¬ 
cluded phenomenologically by adding an ohmic resistor 
and a capacitor in parallel with the junction [resistively 
shunted junction (RSJ) model] Owing to advances in 
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FIG. 2. SNS junction through a normal conductor. We con¬ 
sider two models of the normal part: (i) an arbitrary short 
scatterer with the transmission coefficients {Tc} (Sec. 
and (ii) a quasi-one-dimensional diffusive wire of an arbitrary 
length L (Sec. IIIl. 


nanofabrication technology, current experimental inter¬ 
est has turned towards the study of SNS junctions when 
two super cond ucting terminals are c onnec ted via a nor¬ 
mal regiorP^mUl (including gra i jhen^^nEU^ and the sur¬ 
face of a topological insulatop22H2ll j An important new 
physics in this case is related to the superconducting 
proximity effect,!^ which renders the normal part of the 
junction “partially superconducting”. The strength of 
the proximity effect is characterized by the value of the 
spectral minigap (at x = 0) which can be estimated as 
Eg ~ min(A, li/Tesc), where A is the superconducting 
gap in the terminals and Tesc is the time required for 
an electron in the no rmal region to establish a contact 
with superconductors.l^i^ For a diffusive wire of length 
L with good contacts to superconductors, h/Tesc is of the 
order of the Thouless energy, i?Th = hD/L^, where D is 
the diffusion coefficient. 

Depending on the relation between A and H/tbsc, one 
can distinguish between short (h/Tesc ^ A, with Eg « 
A) and long (h/Tesc ^ A, with Eg A) Josephson junc¬ 
tions. For short junctions, /(x) can be expressed in terms 
of the transmission coefficients of the normal region,!^ 
whereas long junctions require special treatment.^ In 
both cases, the critical current can be written as 

Ic = icGEg/e, (1) 

where G is the normal-state conductance of the junction, 
and Zc ~ 1 is a model-dependent factor. Equation 0 is 
a generalization of the Ambegaokar-Baratoff relation to 
SNS junctions. 

In this paper we make the first step towards the theory 
of the supercurrent decay in Josephson junctions beyond 
the tunnel limit and calculate the life time of the super¬ 
current state due to quantum fluctuations. The normal 
region of the SNS junction will be modeled either by a 
short mesoscopic conductor characterized by an arbitrary 
set of transmission coefficients {Tq,} or by a quasi-one- 
dimensional diffusive wire of an arbitrary length L (see 
Fig. §. The spectral gaps in superconductors are as¬ 
sumed to be equal. We will work in the limit J —>■ Jj, 
(but not too close to Ic so that the WKB approximation 
is still applicable). In this limit free energy E(x) is flat¬ 
tened and X becomes the slowest variable in the system, 
which makes it possible to treat electronic degrees of free¬ 
dom in the adiabatic approximation (see the justification 


in Sec. IVB). 


At low temperatures, T <C Eg, thermal quasiparticles 
responsible for dissipation are frozen out and the junction 
can be described by the imaginary-time effective action 
with the capacitive dynamic term: 


'S'[x('r)] = J 


2(2e)^\dTJ 


( 2 ) 


Such a local-in-time description of the phase dynamics is 
valid only in the supercurrent state with / —> /c. The 
effective capacitance. 


G(x) — C^geom G4x), (3) 

is a sum of the geometric, Cgeom, and intrinsic, C'*(x), 
capacitances of the junction. The latter is determined 
by the response of the Andreev bound states to nonsta¬ 
tionary boundary conditions. These states are known to 
be responsible for carrying the supercurrent .1^ At lowest 
temperatures, their low-frequency dynamics cannot be 
damped (with the kernel |a;| in terms of Matsubara fre¬ 
quency) as it is not related to any dissipation, and thus 
it should be capacitive (with the kernel w^). 

Thus the problem of quantum decay of a nearly crit¬ 
ical supercurrent reduces to calculating the intrinsic ca¬ 
pacitance C* of the junction (evaluated at the critical 
phase X = Xc)- The intrinsic capacitance of the tunnel 
(SIS) junction was obtained in the works of Larkin and 
Ovchinniko’^ and Ambegaokar, Eckern and SchdnP^ 

C^r(x) = (37rG/32A)(l-cosx) (4) 


(here G is the tunnel conductance of the barrier). This 
result has been recently extended to the case of arbi¬ 
trary short junctions by Galaktionov and Zaikin.l^ They 
derived the general expression for G®^°''*(x) in terms of 
transmission coefficients Tq, [see Eq. 0] and found that, 
quite generally, G®^“'*(xc) G/A. We will generalize 
that result further and show that the magnitude of C* (x) 
for an arbitrary SNS junction is determined by the ratio 
of the normal-state conductance G to the minigap Eg in¬ 
duced in the normal part due to the proximity effect [cf. 
Eq. Q for Ic]: 

G.(x)=ac(x)G/Eg. (5) 

The dimensionless capacitance coefficient ac(x) is 
mo del-dependent. At the critical phase difference, 
ac(Xc) is of the order of one, and we will evaluate it 
for an arbitrary short scatterer using Galaktionov-Zaikin 
theory (see Table [l| and for a quasi-one-dimensional dif¬ 
fusive wire with the help of the nonlinear sigma model 
(see Figs. and . 

In the case when the geometric capacitance can be ne¬ 
glected (Ggeom ^ Gsf), we obtain that the instanton ac¬ 
tion, which determines the life time tq oc exp(S') of the 
supercurrent state, depends only on the dimensionless 
conductance of the normal part. 


S 




/x5/4 

Tc) ’ 


( 6 ) 
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TABLE I. Various properties of some model Josephson junc¬ 
tions {Eg = A marks short junctions), and the resulting val¬ 
ues for the capacitance coefficient ac{xc) = Ct,{xc) / {G / Eg). 



Tunnel 

Short wire 

Chaotic QD 

Long wire 

Eo 

A 

A 

A 

3.12 Fxh 

Xc/(7r/2) 

1 

1.255 

1.299 

1.27 

ic 

■kI2 

2.082 

2.186 

3.47 


7r/2 

1.706 

1.791 

2.92 

oc(Xc) 

37r/32 

0.589 

0.689 

0.905 


where Gq = j-Kh is the conductance quantum, and 
the coefficient 7 ^ 1 is model-dependent. Its values for 
several model Josephson junctions are given by 


the expression for the Fano factor in the theory of shot 
noise,™ F = r„(l - T„)/ Ea T’a- ^ 

Application of the scattering matrix approach t o su- 
perconducting hybrid system is a more delicate issue.™^ 
It requires the elements of the normal-state scattering 
matrix to be energy independent on the relevant energy 
scale set by the minigap Fg.™ Since the scattering ma¬ 
trix acquires energy dependence at the scale li/resc, such 
a situation is realized for sufficiently short junctions with 
h/Tesc ^ A and hence Fg « A (a diffusive SNS junction 
with ideal interfaces belongs to this class if Fxh ^ A). 
Then the Josephson current can be found with the help 
of Beenakker formula!^ 

^(x) = ^siny^ ^tanh|A, (9) 


{ 0.874 
1.65 
2.59 
1.83 


(tunnel barrier); 

(short diffusive wire); 

(long diffusive wire); 

(ballistic chaotic quantum dot). 


(7) 


where T is the temperature, and Cq, is the energy of the 
Andreev bound state in the corresponding channel: 

= Ayi-T„sin2(x/2). (10) 


The paper is organized as follows. The zero-tempera¬ 
ture intrinsic capacitance of a short SNS junction with an 
arbitrary normal region specified by its transmission coef¬ 
ficients is discussed in Sec.lIJon the basis of Galaktionov- 
Zaikin formula. In Sec. |III| we present an approach based 
on the nonlinear sigma model which allows us to deter¬ 
mine the intrinsic capacitance of the Josephson junction 
through a diffusive quasi-one-dimensional normal wire of 
an arbitrary length, tracing its behavior from the short 
to long-wire limits. The resulting expression for the life 
time of a slightly subcritical current due to quantum fluc¬ 
tuations is presented in Sec.|IV) where we also discuss the 
limits of applicability of our approach. In the concluding 
Sec. El we summarize our results and discuss them in the 
context of the experimentally observed hysteresis in ex¬ 
tended Josephson junctions. Finally, numerous technical 
details are relegated to several Appendices. 

In our systems of units h = ks = 

II. SHORT MESOSCOPIC CONDUCTOR 
A. Scattering matrix approach 

In this Section we consider the case when the normal 
part can be described by the set of its transmission eigen¬ 
values Ta with the distribution function 

V{T) = [Y.6{T-T^)). ( 8 ) 

OL 

The scattering matrix theory proved to be a powerful 
and intuitively clear method for studying quantum trans¬ 
port in mesoscopic conductors.™ For normal systems, its 
most renowned predictions are the celebrated Landauer 
formula for the dc conductance,™ G = Gq and 


Phase dynamics of a short Josephson junction has been 
recently studied by Galaktionov and Zaikin.™ With the 
help of the Keldysh technique they derived an effective 
action for small phase fluctuations near the equilibrium 
phase X- Slow phase dynamics in the zero-temperature 
limit is governed by the capacitance (§ , where the in¬ 
trinsic capacitance is given by 


- [l-T„sin2(x/2)]”'^' 


2r„ (r„ - 2) sin^ (x/2) 


+ 5 -l- Tq, + 


2 - 2 (l + 2 r^)sin" (x/ 2 ) 
Ta sin^ (x/ 2 ) 


( 11 ) 


Since the resulting expression is sufficiently involved, we 
find it instructive to rederive Eq. (Ill in the Matsubara 


representation. Though the general line of the deriva¬ 
tion is very similar to that of Ref. |44l the absence of an 
additional Keldysh matrix structure makes it possible to 
track the details of the calculation. This procedure sum¬ 
marized in Appendix [A] reproduces Eq. (11) obtained by 
Galaktionov and Zaikin. 


B. Intrinsic capacitance for model junctions 

Equation 0 can be used to evaluate numerically the 
intrinsic capacitance, ^^(x), of various short Josephson 
junctions.™ We calculate it at the critical phase, C'*(xc)j 
for three types of structures, with the superconducting 
terminals coupled through the following links: 

• a tunnel barrier with all Tq, <C 1 , 
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• a short (E’xh ^ A) diffusive wire, with V{T) given 
by Dorokhov distribution, 


V{T) = 


G 


2Gq 

a ballistic chaotic quantum dot,l^ with 


r{T) = 


2G 


ttGq ^T{1 - T) ■ 


( 12 ) 


(13) 


The results along with some characteristics of these 
Josephson junctions are presented in Table|^ We see that 
the capacitance coefficient ac{Xc) defined in Eq. ([^ is 
generally of the order of unity, depending on the partic¬ 
ular distribution of transmission coefficients 'P{T). 


III. QUASI-ONE-DIMENSIONAL WIRE 

Here we calculate the intrinsic capacitance of the SNS 
junction made of a normal diffusive wire coupled to 
superconductors through highly transparent interfaces. 
The scattering-matrix approach used in the previous Sec¬ 
tion cannot be applied to sufficiently long wires (Exh ^ 
A). In this limit, the energy dispersion of the scattering 
matrix at relevant energies E ^ Exh is not negl igible, a nd 
the action does not have a simple form of Eq. (A1 l.l^To 


find C* (x) one then has to use a more general method of 
the nonlinear sigma model and study spatially inhomo¬ 
geneous configurations of the matrix field Q in the wire. 


where prime stands for the derivative with respect to 
r. Assuming perfect NS interfaces, we write the bound¬ 
ary conditions with the antisymmetric choice of super¬ 
conducting phases as Qo je; ± 1/2) = where 

(e) is defined in Eq. (A3). 

In the standard parametrization in terms of the spec¬ 
tral angles 0 e(r) and </ 5 e(r), 

Qo = (''"1 cos ip — T 2 sin (p) sin 9 + cos 6, (16) 

the Usadel equation reduces to two coupled equations: 

(sin^ 9ip'y = 0, (17a) 

6” — 2e sin 9 — {ip')‘^ sin 9 cos 9 = 0, (17b) 

with the boundary conditions at r = ± 1 / 2 : 

0(±l/2) = 9s{e) = arctan A/e, ip{±.l/2) = ±x/ 2 . 

( 18 ) 

The stationary supercurrent is given by 


ttG 


I{x) = -sin^ 6»,(r) (p/(r) 


(19) 


[its conservation is guaranteed by Eq. (17a)]. Tempera¬ 
ture and length dependence of the critical current of dif¬ 
fusive SNS junctions was studied in Ref. Our further 
analysis will be limited to the T = 0 case. 


B. Perturbative expansion near the saddle point 


A. Diffnsive sigma model and Usadel eqnation 

Superconducting proximity effect in a diffusive metal 
can be conveniently d escrib ed by the replica sigma model 
in the imaginary tim d^^ * ^^ ) Its action for the quasi-one- 
dimensional metallic wire of Fig. can be written as 


S[Q] = 


G 


.112 


IhGg y_i/2 


dr Tr 


{VQY - deraQ 


(14) 


where the spacial coordinate r is measured in the units of 
the wire length L, summation goes over Matsubara en¬ 
ergies e = TrT{2n ± 1), and e = e/Exh with the Thouless 
energy Exh = hD/L’^. In general, the field Q subject 
to the constraint = 1 acts as a matrix in the replica 
space, but since the phase difference x carries no replica 
index, the structure of the instanton describing quan¬ 
tum tunneling is trivial in the replica space and we can 
omit the replica index everywhere. Then Q = Qej^^g 2 (r) 
becomes a matrix in the Nambu-Gor’kov space (Pauli 
matrices r^), as in Appendix 

The stationary saddle point QoYi,e 2 ]r) = 27r(5(ei — 
e 2 )Qo(eiu) for the sigma-model action (14) satisfies the 
Usadel equatioii^ 


— {QoQ'o)' + £["^3, Qo] = 0, 


(15) 


In the presence of a time-dependent phase difference 
across the junction, xY)^ the Q matrices in the lead s 
become functions of two time arguments [cf. Eqs. (A4)]: 


Qr(lM,T2) = (n _ ^2)eUxU2)r±2^ 

( 20 ) 

In order to get the action as a functional of x(t) one has 
to integrate out Q(r) in the wire with the boundary con¬ 
ditions (20). Near the criticality, I —)■ Ic, variations of 
x{t) = xyr) — Xc on the instanton trajectory are small 
and hence induced (non-diagonal in energy) deviations 
of Q from the saddle-point Qo can be treated perturba- 
tively. 

Perturbative expansion of the diffusive sigma-model 
near the supercurrent state in Josephson junctions has 
been developed in Ref. 1541 and here we generalized it 
to the case of nonstationary boundary conditions (an al¬ 
ternative approach would be to work on the level of the 
dynamic Usadel equations as it was done in Ref. jSSj). 
In that paper the authors studied the real part of ad¬ 
mittance for the case of the long diffusive SNS junction, 
while the intrinsic capacitance can be extracted from the 
low-frequency expansion of its imaginary part. 

First we present the stationary saddle point 
Qo{r) in the form Qg = U~^tiU with the 
diagonal-in-energy position-dependent matrix 
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U^{r) = Then we 

parametrize small fluctuations near Qq as 

Q = C/-Vi(1 + TT + 1TV2 + ...)^^ (21) 

in terms of the field W satisfying {ti,W} = 0 and = 
— W. The former constraint implies W = cti + dra, and 
we combine c and d into a vector object 


,(r) = KieUO 


2{r). 


( 22 ) 


where the unitary matr^ V is introduced to simplify the 
quadratic action [Eq. (29) below] Its explicit form can 
be found in Appendix |B[ The Pauli matrices acting in 
the (c, d) space will be referred to as .1^ Physically, the 
held w describes soft diffusive modes on top of a super¬ 
conducting state (in the presence of a supercurrent they 
can no longer be classified as difusons and cooperons). 

Expansion of the action in the powers of w reads: 


S[Q] = 5[go] + 


(23) 


where each term is a combination of the bulk and bound¬ 
ary contributions: 

= + (24) 


Since Qo is the saddle point, the linear term is to¬ 
tally due to the boundary, but it will contribute to the 
quadratic-in-a: action [see Eq. (36) below]. The boundary 
contributions are given by 


5, 


5, 


(w^) 

bound 


(w) 

bound 

G I 


iG 


de 

27r 




1/2 

-1/2’ 


4Gq 


dtl dt2 t n 


1/2 

-1/2’ 


(25) 

(26) 


where wl^ = (lUej^ej)!, and the operators b and B read 

b.{r) = , (27) 

, X cos 6e,(pL + cos 9e^(p' , , 

(28) 


c(“") _ 
“^bulk “ 


The bulk contribution is given by 

G rdeide2 
4Gq J L-K ZTT J_i/2 

(29) 

where explicit expression for the effective potential fl can 
be found in Appendix [B| 


f dci de2 


with the boundary conditions (20). It can be written as 
a series in x: 


w = 


■ w 


( 2 ) 


(30) 


with oc x". As we are interested in quadratic-in- 
X terms, it suffices to follow only and Due 

to nonlinearity of the theory, the linear terms may 

influence the equation of motion for the quadratic term 
However, as shown in Appendix]^ this mechanism 
gives no contribution to the effective action. Therefore 
we may consider linear equations of motion obtained by 

Cr(w^) 

varying : 


[-V^ + D,^,ed^)] 


(31) 


and process each order in x independently. 

The boundary conditions for Eq. (311 can be obtained 


by comparing Eq. (201 with the parametrization (21) at 


r = ±1/2. The first two terms are given by 
w'£i?£ 2(±4/2) = 4\L)(D±2)a:ei-6 

= J J?SL)(D± 2 ,e)a:., 

where 

^fl(L)(^i>^ 2 ) = (4=1/2) 





A2) 


ei± 2 , e) — — 
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r.7,/ (± 1 / 2 ) 


(32a) 

(32b) 

(33a) 

(33b) 


and = [6*s(e) ± 9s{e')]/2. In what follows we will 
need only diagonal elements of wG\ so we introduce the 

(O'] (‘2') 

notation ?fl(i)(ei,e2) = Cfl(i)(ei,ei,_e2). 

With the boundary conditions 
Eq. (31) can be written in the form: 


the solution of 


-,(i) 






e2'^€2 — €l 5 


(34a) 

(34b) 


where the functions ^ei^es ("c) and Cei?e 2 (’") obey the same 
Eq. (31) as and satisfy the boundary conditions: 


/!;,>„(±l/ 2 ) = 4 ',/,(e„e,), 

{!±(±l/ 2 ) = 4h)('i.'2)- 


(35a) 

(35b) 


To get the quadratic-in-x part of the action, we substitute 
Eqs. (34) to Eq. (23): 


C. Effective capacitance evaluation 

The effective action for the phase variable x(t) is de¬ 
termined by the saddle-point trajectory of the action (23) 


S'(2) = ±S'(“)[u;(2)]^ 


(36) 


which takes the form of Eq. ( |A5[ ) with the kernel /C(e, w) 
after the substitution e = (ei ± e 2)/2 and ui = ei — £ 2 - 
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Then the intrinsic capacitance should be determined 
from Eq. (A6). To find it numerically, we expand the 

matrix il.e+uj/ 2 ,e-uj/ 2 {r) and the solutions Cf+^/ 2 ,£-a;/ 2 (^) 
in powers of ui as 


^e+u,/ 2 ,e-uj/ 2 (r) = ^ fi/c (e, r), 




cl/2, e-a;/2W = H C 




(37) 

(38) 


and write down the system of linear equations for the 
first three coefhcients: 



A/i?Th 


(_v 2 + = 0 , 

(_V2 + (39) 

(_V2 + - r!2e^"\ 


which should be supplemented by the boundary condi¬ 
tions obtained from the expansion of + w/2,e — 

uj/2) in uj. Extracting jdui^ from Eq. (36), we 

obtain 


i9^/C(e, w) 


9^2 


= > 




1/2 

- 1/2 


+2i:fr"v{<b 


k=0 


1/2 

- 1/2 

(40) 


[the term originating fro m th e S 2 part of Eq. (28) equals 
zero]. According to Eq. ( |A6| ) , the intrinsic capacity of the 
junction is determined by the integral of Eq. (40) over e. 


D. Notes on numerical evaluation 


The first step of numerical simulation is to find the so¬ 
lutions for the Usadel equations ( |T7| . One way is to use 
the explicit expression in terms of elliptical functions,!^ 
that reduces to finding two constants from a system of 
algebraic equations numerically. However in order to ob¬ 
tain the functions O^ir) and (/>e(r) numerically we find 
it more convenient to follow the other way and to simu¬ 
late the original differential Eqs. 0 directly. Obtained 
spectral angles are used to determine the matrices Og, 
ill, and ^2 in Eq. (37). Then the system (39) is solved 
numerically, the solutions are substituted to Eq. (40) and 
the integral over e in Eq. (A6) is calculated with a proper 
energy grid. 

As a check of the numerical method, it is instructive 
to calculate / deJC{e,0). It is easy to see that it should 
be proportional to the derivative of the current-phase 
relation dijdx- By checking that this quantity indeed 
crosses zero at the critical phase Xc we can provide an 
independent test for our numerical calculations. 


FIG. 3. The capacitance coefficient at the critical phase, 
oc(Xc) = C*{Xc)/{G/Eg), of the Josephson junction made of 
a finite-size wire vs. A/Axh- It is changing from a^°'^*(Xc) = 
0.589 for short wires to ac“®(Xc) = 0.905 for long wires. 



FIG. 4. The capacitance coefficient ac(x) = C*{x)/{G/Eg) 
of the long (Axh A) Josephson junction through a quan¬ 
tum wire vs. the phase difference across the junction, y. The 
critical phase, Xc = l-27(7r/2), is marked b y th e arrow. The 
dashed line shows the asymptotic behavior (41) near y = tt. 


E. Results 

The resulting dependence of the capacitance coeffi¬ 
cient Qfc(Xc) = C^{xc)/[G/Eg) evaluated at the critical 
phase difference, Xc, on the ratio A/Axh cx is pre¬ 
sented in Fig. (for the dependence of Eg on A/Exh 
see, e. g., Fig. 5 of Ref. |54]). In general, aciXc) is of 
the order of unity, varying from ap°‘'*(xc) = 0.589 for 
short wires to Q;p”®(xc) = 0.905 for long wires [in that 
limit, Eg = 3.12 Exh (Ref. [571) and h = 10.82 GExh/e 
(Ref. |24|)]- The maximal value of a^^^(Xc) = 0.958 is 
achieved at A/Exh = 5.8, corresponding to the wire 
length L = 2A^~D]~K. 

For completeness, in Fig. we present the phase de¬ 
pendence of the intrinsic capacitance coefficient ac(x) = 
G*{x)/{G/Eg) in the long-wire limit. In the vicinity of 
X = TT, where the gap closes, the intrinsic capacitance 
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diverges as 


c^r^ix) 


0.55 Gn 
(tt - x)2 Eg ■ 


(41) 


A similar behavior obtained for short diffusive junctions 
in Ref. |44] is attributed to the presence of almost open 
channels (Tq —>■ 1) in Dorokhov distribution (12). 


IV. 


LIFE TIME OF THE SUPERCURRENT 
STATE 


If the geometric capacitance is sufficiently small, 
Cgeom ^ C'*(Xc), the instanton action ([46| ) can be sim¬ 
plified. Substituting C^^Xc) from Eq. (f^ we arrive at 
Eq. ([^, with the model-dependent parameter 


7 = 


12x23/4 ^5/4^ 


5tt 


(,v)3/4 


(47) 


For some model cases the coefficient 7 is given by Eq. 0- 


B. Applicability of the theory 


A. Instanton action 

Decay of the dissipationless supercurrent state gov¬ 
erned by the action ([^ is equivalent to a quantum- 
mechanical tunneling of a massive particle under a bar¬ 
rier. Near the criticality, at / —> Ic, the potential barrier 
AF{x) = F{xc + ic) — F{xc) can be approximated by a 
cubic parabola: 






*c( t - — )a: - —X 


(42) 


where ic is defined in Eq. (§,and*" = -/"(xe)/(GE,/e) 
is the dimensionless curvature of the current-phase char¬ 
acteristic at Xc- Parameters ic and i" for some model 
Josephson junctions are listed in Table |l] The instanton 
trajectory for the potential (42) takes the form 


x{t) -Xi = 


cosh^ {ijjpT/2) 


, 1/2 


(43) 


where LOp is the plasma frequency describing small phase 
oscillations near the minimum of AF{x). In the short- 
wire limit, ujp was calculated in Ref. 1441 and the general 
expression is given by 


UJp = 23/4 (/^i") 1/4 


I EgG 
G(Xc) 


(>4) 


1/4 


(44) 


The life time of the subcritical current state in the 
WKB approximation can be estimated as: 


To « {2Tr I ujp)e^ 

where S is the instanton action, 


12x23/4 -5/4 ^C{Xc)EgG 

57r (i")^/^ Gq 


S = 




5/4 


(45) 


(46) 


This expression generalizes the result for the tunnel limitP 
to the case of an arbitrary Josephson junction. 

Equations ([44|-(46) supplemented by the expressions 
^ and ([^ provide full description of the supercurrent 
decay due to quantum fluctuations. The numerical coeffi¬ 
cients ic, i'c and adXc) for a number of model Josephson 
junctions are summarized in Table [ij 


In our analysis we rely on the adiabatic approxima¬ 
tion for xd) justified at J —>■ Ic- In this limit, variation 
of the phase is weak and small, that allows to expand 
the action. On the other hand, in the very vicinity of 
Ic, the instanton action 5” < 1 and the WKB method 
fails. Therefore the present theory is applicable as long 
as [GQ/G(xc)AgG]^/3 ^ 1 — //A ^ 1. If the geometric 
capacitance is negligible that reduces to 


(Gq/G)4/3 < 1 - ///, < 1. 


(48) 


Finally, we discuss the condition on the temperature 
range when the T — 0 description of the phase tunneling 
is applicable. The principal limitation is related to the 
crossover to thermal decay with the activation exponent 
where SF ~ {G/GQ)Eg{l- 1 /Ic)^G ig the height 
of the free energy barrier. Quantum description is appli¬ 
cable as long as SF/T 3> S, which translates into the 
constraint [C{xc)T'^/GEg]"^ <C 1 — I/Ic- In the impor¬ 
tant limit of a small geometric capacitance, this condition 
reduces to a simple inequality 


{T/Eg)^ « 1 - l/lc- 


(49) 


The condition ( |49[ ) also guarantees that the phase action 
can still be written in a capacitive form 0 neglecting 
quasiparticle damping effects. Indeed, dissipative terms 
in the action can be roughly described by a large shunt¬ 
ing resistance R{T) ~ , which has no effect on 

phase tunneling since the corresponding RG time is much 
larger than the period of plasma oscillations. 

To estimate the predictions of our theory consider 
a 1 /rm-long SNS junction similar to that fabricated in 
Refs. E51 E5l and E51 With the parameters G“4 ~ 5S1, 
Erh ^ 30mK, and A/Exh ~ 10^, the intrinsic ca¬ 
pacitance is expected to be of the order of 10 pF. The 
condition (48) determines a wide range of allowed bias 
currents: 10“3 1 _ I/I^ <c 1 , in which the pre¬ 

dicted life time (when the intrinsic capacitance domi¬ 
nates) varies from nanoseconds to practically unlimited 
values. For example, the decay time tq ^ 1 s is achi eved 
at l — I/Ic ~ 10“^ . According to the constraint (49), the 


quantum tunneling regime is realized in the sub-100 mK 
temperature range. 

Our calculation was performed in the assumption of 
ideally transparent SN interfaces, whereas experimen¬ 
tally studied junctions may have a finite conductance Gt 
















of interfaces. Our theory is valid as long as Gt ^ G, but 
it can be straightforwardly extended beyon d this limit 
by adding the corresponding boundary term^^^lMI to the 
sigma-model action (14|. In the opposite limit, Gt ^ G, 
the normal region can be treated as a chaotic quantum 
dot as it was done in Ref. EH 


V. DISCUSSION AND CONCLUSION 

In the present paper we have found the life time of a 
slightly subcritical dissipationless supercurrent state in 
an SNS Josephson junction due to quantum fluctuations. 
At low temperatures, the superconducting phase differ¬ 
ence across the junction behaves as a massive quantum- 
mechanical variable, with the mass determined by the 
sum of the geometric (Ggeom) and intrinsic (G*) capaci¬ 
tances of the junction. While the former describes charg¬ 
ing effects, the latter is associated with the dynamics 
of the Andreev bound states. We obtain that, generi- 
cally, the intrinsic capacitance is expressed through the 
normal-state conductance G and the spectral minigap Eg 
induced in the normal part of the junction by the relation 
C*{x) = OLc{x){G/Eg), where ac{x) -- 1 is a model- 
dependent coefficient. An analogous expression for G* 
in the long-wire limit was conjectured in Ref. 1251 based 
on the requirement that the i?G* time of the junction is 
governed by h/Eg. 

We determine the model-dependent numerical factor 
ac(Xc) for two classes of Josephson junctions: an arbi¬ 
trary short scatterer described by the set of its transmis¬ 
sion eigenvalues (see Tablejl]) and a diffusive metallic wire 
of arbitrary length (see Fig. [^. In the former case we 
rely on Galaktionov-Zaikin expression ( |11[ ), rederived in 
Appendix]^ In the latter case we employ the nonlinear 
sigma model formalism. Depending on a particular junc¬ 
tion, the capacitance coefficient aciXc) can vary by a 
factor of three that may have a pronounced effect on the 
quantum tunneling rate due to its exponential sensitiv¬ 
ity to the junction parameters. For long junctions with 
negligible geometric capacitance, the decay rate is deter¬ 
mined only by G and proximity to criticality, as given 
by Eq. (§. Measuring the coefficient 7 in that relation 
and comparing it with the results 0 may provide an in¬ 
dependent tool to determine the type of the Josephson 
junction. 

Our analysis of the quantum decay of the supercurrent 
is limited to low temperatures, T Eg. For higher tem¬ 
peratures, finite population of quasiparticle states leads 
to dissipative phase dynamics, and the problem should 
be treated in the spirit of Ref. [3l with an additional com¬ 
plications due to inapplicability of the simple RSJ model, 
as well as possible issues on inelastic relaxation and ther- 
malization in the normal part of the junction.l^ 

Capacitive phase dynamics is realized only in the 
nondissipative supercurrent state near the criticality. 
Once it is destroyed by fluctuations, the junction switches 
to the resistive branch and the phase dynamics becomes 


dissipative (like for high temperatures). Though we can¬ 
not access this regime, our results can be used for a quali¬ 
tative description of the hysteretic behavior observed ex¬ 
perimentally in lateral junctions.l^^ll^ Assuming that an 
oversimpified RSJ model can qualitatively describe the 
SNS junction with the choice of the resistance R ~ G~^ 
and capacitance G ^ G^, one finds that the McCumber 
parameter Pc = 2eIcGR^ (which is the square of the 
quality factor Q = cOpRC) is generally of the order of one 
[formal substitution R = G~^ and G = G*(xc) yields 
Pc = 2 zcac(xc)) which is larger than 6 for long wires]. 
In this model, large McCumber parameter, /Jc > 1, is 
required for a hysteretic behavior. Though for long junc¬ 
tions the geometric capacitance is too small to explain 
the hysteresis, an account of the intrinsic capacitance 
provides sufficiently large Pc ^ 1 necessary for observing 
a hysteretic behavior.l^ Recently it was demonstrated 
that low retrapping current is a consequence of elec¬ 
tron overheating in the normal region, when the electron 
temperature can be several times larger than the bath 
temperature.^ The retrapping current is then identified 
with the critical current at the elevated electron temper¬ 
ature, explaining large hysteresis. By contrast, we would 
like to emphasize that the intrinsic junction capacitance 
due to dynamics of the Andreev bound states may itself 
lead to a hysteretic behavior, even for a perfect thermal 
contact with the environment. 
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Appendix A: Intrinsic capacitance of a short 
Josephson junction 


In this Appendix we rederive Calaktionov-Zaikin for¬ 
mula (11) in the Matsubara formalism. 


1. Matsubara action 


An arbitrary scatterer with an energy-independent 
scattering matrix (short-wire limit) sandwic hed betw een 
two terminals can be described by the actiorPSEHini 



^Trln 


1 


4 


{Ql - Qr) 



(Al) 


where Ql and Qr are the quasiclassical Creen functions 
in the leads. Q[ei,t 2 ) is a function of two energy ar¬ 
guments and acts as a matrix in Nambu-Cor’kov space 
(N). When multiplying Q-matrices and taking trace, in¬ 
tegration over ei and 62 is assumed to be done. The time 
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representation is defined in the conventional way as 

Qiri,r,)= |^Q(ei,e 2 )e—(A2 ) 

In a stationary uniform superconductor, Qs depends 
only on the time difference, QsiTi,T 2 ) = (5 s(ti — T 2 ), 
and hence is diagonal in the energy representation, 
Qs(ei, £ 2 ) = 27r(5(ei - e 2 )( 5 s(ei), where 




1 ( e 

Ve^ + A 2 "Wn’ 


(A3) 


and If is superconducting phase. 

Now we proceed to the derivation of the capacitive 
term in the action ([^. We assume zero superconducting 
phase on the left lead and a time-dependent phase x(t) 
on the right lead (see Fig. [^. We assume that x(r) = 
Xi + x{t), where x(t) is a small and slow function of 
Matsubara time r: 

Ql(ti,T 2 ) = Qs^in - T 2 ), (A4a) 

QH(n,T2) = (A4b) 


where is the Pauli matrix in the N space. Substituting 
these expressions into Eq. ( |A1[ ) we obtain the action for 
the phase difference x(t). 

In the limit / —>■ /c, the free energy barrier protect¬ 
ing the supercurrent state X = Xi is small, and one can 


expand the action in powers of x(t) = xi^) — Xc- The 
quadratic term can be written in the form 


S(2) 


G 

4^ 


doj 

27T 


^ XujX-^ lC{e,u}), 

ZTT 


(AS) 


with the kernel A(e,w) calculaed below. Also, the same 
limit guarantees that phase dynamics is slow and thus 
can be described (at T = 0) by the term x^ in Eq. (|^, 
with the intrinsic capacitance of the junction given by 


C* = ttG 


/ 


de d'^JC{e,uj) 
27r diu'^ 


(jj—0 


(A6) 


2. Derivation of the kernel K, 


Here we calc ulate the kernel IC{e,uj) starting with 
the action (Al) for an arbitrary distribution of Ta (for 
brevity we omit brackets that denote averaging over Ta). 
We write Qr as 

QRiTi,T2) = e“("^A^AQfl(ri,r2)e-“(-2AA2^ (A7) 

where Qr{ti,T 2 ) = — T 2 ) corresponds to the 

stationary phase xi- Then we insert expressions (A4) 
into Eq. (AI) and expand the action in powers of x{t) = 
Xix) - Xi- 


The first variation of the action with respect to x reads: 


= — Tr 


AQ 


8 1 - {Ta/^)AQ^ 


TzX, 


.A1 = - EE T 


Cf k—0 


Ta 


fe + 1 


Tr AQ 


2fe+lr 


zXj Q r\ , 


(A8) 


where we denote AQ = Ql — Qr and {x)^^, = x{t)5tt'■ 

In the following expression for the second variation we will omit for simplicity terms that contain two x operators 
without any Q-matrix in between, since such terms will give no contribution to the intrinsic capacitance. Also after 
taking the variation, we put x{t) = 0 as we are interested in the second variation near the stationary solution. The 
result can be written as 

5'^S = Tr< ^ AQ\tzX, Qr]AQ'^^~'-TzxQb. + AQ^’^TzxQrTzxQr T AQ'^^+^TzxQrTzX 

oi k^Q ^ ^ I 1^0 

2k 

- ^ AQ^ItzX, Qr]AQ‘^^~^QjiTzX + AQ‘^^~^^TzxQrTzX + AQ‘^^QrTzxQrTzX 
1^1 



where AQ = — Qr is diagonal in the energy representation: 


2 Asin(x/2) 
Ve 2 + A 2 


AQ(e) 


q = Tz: sin(x/2) -f Ty cos(x/2). 


(AlO) 
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Switching to the energy representation, we arrive at the action (A5) with the kernel 


/C(e,a;) =) I ^ 2X2fc,i [trT^Qfi(e)T^Qfl(e') - tr 1] 


fc =0 


2k 


. Z=0,2,., 


2k-l 


+ X2k,i [‘^^'^r:,QR{e)qT^Qii{e')q-ivT:,Qii{e)qQR{t)T^q-tvT^QR{e')qQii{e')T:,q] 

+ -^2fc,0 + ^2k,2k + -A2fc+1,0 tl' QTzQR{ f)Tz + A2fc+i^2fc+l tr QTzQR{e')Tz > , (All) 


where e' = e — ui, and 


Evaluating traces one gets: 


Xnl = 


[2Asin(x/2)]^ 


(e2 +A2)'/2(e'2 + 


rp py OO / rp \ ^ { 

\- Y. 

a k—0 ^ ' 


VA2 + eVA2 + e'2 


2k-l 

Y. 2X2k,i ( 2 


/=0,2,... 

ee' — A^ cos x + A^ cos x cos x 


(A12) 


v/A2 + e2v/A2 + e'2 e2 + A2 e'2 + A2 

+ 2A2fc,0 + ‘2X2k,2k + A2fc+2,1 + X2k+2,2k+l )■ • 


(A13) 


Performing summation over I and k one obtains for the kernel X{e,uj): 


= Y 


TaGg 


A2 


+ 2 


8 G e2 + A2 (1 - 74 sin^ f) 

-v/A2 + + e'2 f ^ gg' 


-4- 


+ A^ 


1 + 


£2 _ g/2 VA2 + e2VA2 + e'2 


- 2 - 


(2A sin 


- A2 


cos X 


+ A^ cos X + A^ cos x 


VA2 + eWA^ + e'2 

+ {e O e'} . (A14) 


v/A2 + e2v/A2 + e'2 e2 + A2 e'2 + A2 

Finally, expanding this expression to the second order in uj and integrating over e we arrive at Eq. 0 obtained by 

Galaktinnnv and Zaikinl^ 


Appendix B: Explicit form of the matrices V and fl we omit the spacial coordinate r for brevity)!^ 

^ei,e2 “ *^ei,e2 Pei,£2 A 7e2)^3 

+ Pei ,£2 sin(77<,, + 77^2)Ei, (B3) 

where 


The rotation by the matrix V in Eq. (221 is introduced 
in order to get rid of the first derivative in the bulk action 

•^buik • It® matrix elements are given bjEl 

= e*S2[C.iM+C.2M]/2^ (Bl) 

where are Pauli matrices in the (c, d) space, and 


Ce(r)=-/ ds cos 9 eis)ipeisy. 

JO 


(B2) 


( ^4 

After such a rotation, acquires the form (29), with 
the matrix potential f2e^^e2(r) given by (here and below 


aei,e 2 = £1 cos Pei - + (sin P^^)^] /4 

+ £2 cos Pe2 - [P(^ + (sin 9^^tpyy] /4, (B4) 

and 

Pei, £2 = + (sin^'eiP'ei)^y^fi '(2 + (sin ^£2 P'e 2 

(B5) 

and the function ri^{r) is given by (as usual, e = e/Exh) 

sin^ Pe(s) I75'e(s) 


77 e(r) = —2e / ds 


0 P'2(s) + [sinPe(s):p((s)]2 


(B6) 
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Appendix C: Irrelevance of nonlinear terms in the 
action (231 


In this Appendix we show that for finding the capac¬ 
itive effective action (pi) it is sufficient to use the linear 

equation (311 [resulting from varying only and ne¬ 

glect higher order terms in Eq. (23). 


The only possible extra contribution to the quadratic- 
in-a: part of the action can originate from the influence 
of on the equation of motion for arising from 
cubic nonlinearity in the action (231. That w ould induce 
a correction to which satisfies Eq. (31) with a 


nonzero right-hand side: 






( 2 ) 

ei€2 




* foifc f(l) r(l) 
27r 


(Cl) 


where the indices i, j, and k refer to the (c, d) space, and 
the operator h may contain derivatives with respect to 
r. The crucial point is that the boundary conditions for 
are trivial: 


5f(2)(±l/2) = 0, 


(C2) 


since they are already satisfied by [see Eq. (35b)]. 
When substituted into could have produced 

a correction to K.. But since [Eq. ([25 1] is propor¬ 
tional to w at the boundary, condition (C21 guarantees 
the absence of corrections to 1C. 
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